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ABSTRACT 

The Lya absorption systems observed in the spectra of QSOs are hkely to 
possess bulk pecuhar velocities. The free electrons in these systems scatter 
the microwave background and distort its spectrum through the kinematic 
Sunyaev-Zel'dovich effect. I calculate the temperature fluctuations of the 
microwave sky due to variations in the number of Lya systems along different 
lines-of-sight throughout the universe. The known population of absorbers out 
to z ~ 5 introduces anisotropies on angular scales ^ 1' with an rms amplitude of 
order AT/T ^ 10~^(f2LyQ/0.05)(t'4QQ)"'^/^, where ^Lya is the cosmological density 
parameter of ionized gas in Lya absorption systems, and {vIqq}^^'^ is the rms 
line-of-sight peculiar velocity of these systems at ^ ~ 3 in units of 400 km s~^. 
Detection of this signal will provide valuable information about the cosmic 
velocity field and the gas content of Lya absorption systems at high redshifts. 



Subject headings: cosmology:theory-quasars:absorption lines-cosmic microwave 
background 
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1. Introduction 

Sunyaev and Zel'dovich showed that Compton scattering off free electrons would 
distort the spectrum of the cosmic background radiation (CBR) either if the electrons have 
a different temperature than the radiation or if they possess a bulk velocity relative to the 
cosmic frame (see Sunyaev & Zel'dovich 1980 for a review). By now, the thermal effect has 
been seen in many clusters of galaxies (Birkinshaw 1993; Rephaeli 1995), and observations 
to measure the kinematic effect are underway (see discussion in Haehnelt & Tegmark 
1995). The integrated level of distortion introduced to the CBR spectrum by intergalactic 
gas is currently unknown, although it may be considerable in popular cold dark matter 
cosmologies (Colafrancesco et al. 1994, 1995; Persi et al. 1995). The perfect blackbody 
spectrum detected by the COBE sateUite sets a tight hmit on the total Comptonization 
parameter of the universe, < 1.5 x 10~^ (Fixsen et al. 1996). Although this limit can be 
used to rule-out exotic cosmological scenarios (Levin, Frccsc, & Spcrgel 1992; Wright et al. 
1994), it still allows for a significant contribution to the observed CBR anisotropics by the 
Sunyaev-Zel'dovich (SZ) effect. As future satellite experiments are currently being designed 
with the hope to map the microwave sky on angular scales ;^ 10' (Janssen et al. 1996), it is 
important to characterize the unavoidable contribution from intergalactic scattering to the 
expected anisotropics. 

The intergalactic medium at redshifts 2; ^ 5 is highly ionized and inhomogcncous. 
Direct evidence for condensations of gas comes from the absorption lines observed in QSO 
spectra (see reviews in Blades, Turnshek, & Norman 1988). Close QSO pairs and multiple 
images of gravitationally lensed QSOs indicate that the transverse extent of these absorption 
systems is ~ 0.1-lMpc (Bechtold et al. 1994; Dinshaw et al. 1994, 1995; Fang et al. 
1995; Smette et al. 1995). The neutral hydrogen density in these systems is thought to be 
controlled through photoionization equilibrium by an ionizing radiation background. The 
existence of this background radiation has been established by the "proximity effect" , with 
the background flux estimated from the size of the ionization region surrounding bright 
QSOs (e.g., Bajtlik, Duncan, & Ostriker 1988; Lu, Wolfe, & Turnshek 1991; Bechtold 1994). 
Hydrodynamic simulations of popular cosmological models tend to recover the essential 
properties of the observed absorption systems and imply that these systems are a generic 
consequence of early structure formation in the universe (Cen et al. 1994; Zhang et al. 
1995; Hernquist et al. 1996; Miralda-Escude et al. 1995). 

Like other objects in the universe, the Lya absorption systems are expected to develop 
bulk peculiar velocities in response to the gravitational field in their vicinity. Thus, as the 
free electrons in an absorption system scatter the microwave background, they would distort 
its spectrum through the kinematic Sunyaev-Zel'dovich effect. For a single absorption 
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systems with a line-of-sight peculiar velocity of \v\ ^ 100 km and a characteristic 
temperature in the range T ~ 10^~^ K, the kinematic SZ effect is much larger than the 
thermal SZ effect since |f/c| ^ 10 x (2/cT/meC^), where rrie is the electron rest mass. This 
is the opposite situation from that in clusters of galaxies. The Rayleigh- Jeans distortion 
is then (AT/T) ^ —rv/c, where r is the optical depth of the Lja system to Thomson 
scattering. The contributions from positive and negative velocities of individual systems 
would tend to cancel out on average, and the integrated CBR spectrum over the whole sky 
would therefore show little evidence for this distortion. However, in small regions of the 
sky, of order the angular scale occupied by a single Lya absorption system 1'), there 
could be large rms amplitude of fluctuations due to incomplete cancelations between a small 
number of absorption systems. In this Letter, I calculate these fluctuations and show that 
their detection could provide important new information about Lya absorption systems at 
high redshifts. 



2. Calculation 

Our discussion is based on the assumption that the intcrgalactic gas at high-redshifts 
is in photoionization equilibrium with a UV radiation background. The time required 
to establish photoionization equilibrium, ~ 10'^J2i^ yr, is much shorter than the Hubble 
time, for the Lyman-limit intensity inferred from the proximity effect (Bechtold 1994), 
Jf = J21 X 10~^^erg cm~^ s"^ }iz~^ sr~^. We focus our attention on systems with column 
densities below the Lyman hmit, < A^max = 10^^ cm"^, because thicker systems are 
opaque to photoionizing photons and are therefore composed primarily of neutral hydrogen 
that does not interact with the microwave background. 

Let us flrst flnd the rms amphtude of fluctuations in the optical depth of the universe 
due to the dumpiness of the intcrgalactic gas. Most of this gas is known to be ionized 
and its neutral hydrogen component serves as a tracer of the fluctuations in the electron 
density. At photoionization equilibrium, the dumpiness in the electron density is directly 
proportional to the local neutral hydrogen density, oc nn- Thus, overdense regions show 
up as absorption lines in QSO spectra and their statistical properties can be inferred from 
observations. We define Ug to be the mean electron density within an absorption system of 
thickness £ and HI column density A^. At photoionization equilibrium, 

nl£ = TjN (1) 

where 77 30J2i(r/5 x 10^ K)°-^ cm-=^ (Spitzer 1978; Peebles 1993). For the redshift 
interval of interest 2 ^ 2; ^ 5 we assume that J21 and therefore 77 are constant in time (see 
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discussions in Bajtlik et al. 1988 and Bechtold 1994). The optical depth of an absorption 
system to Thomson scattering can then be obtained from equation 

r{N) = arnJ = ar^vNlf'^ (2) 

where is the Thomson cross-section, and ^(A^) is the ensemble averaged thickness of all 
absorption systems with column density N . For simplicity, we parametrize I = BN'"^. 
While low-column density systems with N ~ 10^^ cm^^ may have a line-of-sight thickness 
I ^ 10^ kpc (Ranch & Haehnelt 1995; Bechtold et al. 1994; Dinshaw et al. 1994, 1995; Fang 
et al. 1995; Smette et al. 1995), damped Lya absobers with ~ 10^^ cm~^ are probably 
proto-galactic disks with a characteristic thickness £ ~ 10°~^ kpc. This modest variation in 
thickness over many decades in column density implies a relatively small value of a ~ 0.2. 
As shown later, the value of B does not enter into our final result. 

The probability of intersecting a Lya absorption system with an HI column density 
between A^ and A^ + dN in a redshift interval dz around 2 ~ 3 is well fitted by a power-law 
(Carswell et al. 1984; Tytler 1987), 

dP \ 



dNdz) , 



= AN->^, (3) 

=3 

where A = lO^-^^^-^ cm^/^+a] and p = 1.51 ± 0.02 (Sargent, Steidel, & Boksenberg 1989). 
Most of the fluctuations in the optical depth of the universe are contributed by systems 
close to the Lyman-limit. For these systems we extrapolate equation (^ to other redshifts 
with the scaling, 

where 7 = 1.5 ± 0.4 (Stengler-Larrea et al. 1995; see also Storrie-Lombardi et al. 1994). 
Since the systems of interest are separated by large redshift intervals along a given line of 
sight (e.g., there is typically one Lyman-limit system per unit redshift), we can safely ignore 
any correlations between them (see also Sargent et al. 1980; Webb & Barcons 1991) and 
assume that they are randomly distributed in redshift and column density according to the 
probability distribution in equations (|^) and @. 

We can now calculate the rms level of fluctuations in the optical depth out to z = 5 due 
to the finite number of absorption systems along random lines of sight. In each infinitesimal 
(AA^, Az) bin there is either one or no Lya absorption system. The probability of having a 
system with an optical depth t{N), is p = {dP/dNdz)ANAz <^ 1. Therefore, the variance 
of T in each bin is pr"^ — {pr)^ ~ pr"^. The integrated variance is then obtained by summing 
the independent contributions of all infinitesimal bins, 

'AAr]ABaf 



f5 /-Afmax dP 

(Ar') = dz dN^:— X 



2-(3-a 



{N'^.t^ - AT^Tf-") , (5) 
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where the second equahty follows from equations (^-(||) with 7 = 1.5. As expected, the 
total rms amplitude of fluctuations grows in proportion to the square-root of the number of 
Lya systems, namely (Ar^)^/^ oc A^l'^ . The constant coefficient in square brackets depends 
on a variety of uncertain parameters, such as J21, T, or B. We therefore prefer to trade 
these parameters in favor of a single quantity, the density parameter of ionized gas, whose 
value can be related to Big-Bang nucleosynthesis. We achieve this conversion by calculating 
the derivative of the mean optical depth of the universe from equations (H) and @, 



dz 



2=3 



,.r ( dP \ 

dN — — X r 
\dNdzl^^^ 



A{r]By/^aT 
3/2- p- a/2 



\ max 



■P-a/2 



n: 



3/2-f3-a/2 



(6) 



Since the mean optical depth is invariant to dumpiness for an optically-thin medium, we 
can also express d{T)/dz in terms of the density parameter of ionized gas in absorption 
systems below the Lyman-limit, ^Lya, 



dz 



cdt \ 

^''^^''^irz) 



5.5 X 10" 



F(0.3,0.7,3) 



0.05 



(7) 



where /iso is the Hubble constant in units of 50 km s ^ Mpc ^, i7 and are 
the density parameters of matter and the cosmological constant respectively, and 

F(Sl, ^A, Z) = {\^ ^)V[^(1 + + (1 - - ^^a)(1 + zf + fi; 
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The density of baryons predicted by Big-Bang nucleosynthesis could either be in the 
form of ionized gas, neutral gas, or compact objects such as stars. The neutral hydrogen 
can be weighted by the frequency of absorption systems in the spectra of QSOs; this 
yields JIhi ~ (1 — 4) x 10~^ at a redshift 2; ~ 3, a value which is comparable to the mass 
density of luminous stars in the local universe (Lanzetta, Wolfe, & Turnshek 1995). The 
metalicity of the densest absorption systems which are thought to be the progenitors of 
present-day disk galaxies (the so-called damped Lya systems), is only a tenth of the solar 
value (Lanzetta et al. 1995; Pettini, Lipman, & Hunstead 1995; Sembach et al. 1995), 
implying that star formation is only at its infancy at redshifts z ~ 2-5. If we therefore 
assume that the stellar mass fraction is small with a density parameter fij, ^ then 
most of the baryons in the universe may be in the form of ionized gas at 2; ~ 2-5. Based on 
the baryonic density parameter predicted by standard Big-Bang nucleosynthesis (Walker 
et al. 1991), f2b 0.05/i^o^, it is then plausible to substitute ^Lya ~ 0.05 in equation (^. 
Indeed, hydrodynamic simulations indicate that more than 85% of the baryons reside in 
absorption systems below the Lyman-limit (see Fig. 23 in Miralda-Escude et al. 1995). 



By multiplying the square-root of equation (|^) with the ratio of equations (|^) and 
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we get 

nLyaho \ (-3/2 + f3 + a/2) [iV£-/^-") - iVa;"-"^]^/^ 
0.05 ; [A{2 -(3- a)] V2 [iv(„3/2-/.../2) _ ^m-^-^m^ ' 

(8) 

For A = lO^-i cm^ p = 1.5, ^]Lya/i50 = 0.05, n = 0.3, = 0.7, a = 0.2, A^^ax = 10^^ cm~2, 
and A^min = 10^^ cm~^, we get (Ar^)^/^ = 7 x 10~^. This corresponds to an rms amplitude 
of fluctuations of order 3% in the mean optical depth of the universe out to 2; = 5. This 
result has a weak dependence on uncertainties in the values of most of the above parameters 
(to a first approximation, (Ar2)i/2 qc A-^/^jyai^jyaiS)^ 

We finally wish to calculate the resulting fluctuations in the microwave background 
temperature due to the kinematic Sunyaev-Zel'dovich effect of the Lya forest. For this 
purpose, we assume that the optical depth r and the peculiar velocity v of an absorption 
system are uncorrelated and then repeat the calculation that led to equation (|). Note that 
we define v to be the net bulk velocity of an absorption system, and so we average over 
internal velocities within the system. The resulting rms amplitude of fluctuations in the 
microwave temperature is, 

(^) - « 10- (^) P) 

where (f ^) is the dispersion in the line-of-sight peculiar velocities of absorption systems, 
and the second equality was obtained for the set of parameter values chosen after equation 
(§). The dispersion (w^) is averaged over redshift with a weight factor of (1 + . The 
anisotropy signal expressed in equation (|^) varies on the angular size of a Lya absorption 
system at 2; ~ 5, 6'Lya ~ 1' x {i±/500h^Q kpc), and dechnes on larger scales roughly as 

In the local universe, the peculiar velocity field of galaxies shows a characteristic 
one-dimensional dispersion of order ~ 500 km s~^ (Strauss, Cen, & Ostriker 1993) and 
a coherence length of order a few tens of Mpc (Gorski et al. 1989; Strauss & Willick 
1995). Both linear theory (yielding f oc (1 + z)~^'^ for Q = 1 and a more moderate 
decline for i7 < 1) and nonlinear dynamics predict that the velocity dispersion would 
scale down only by a modest factor of order 2 for 2; ~ 3. It is therefore plausible to 
consider (w^)^/^ ~ 300 km s~^ in equation (|^). The redshift evolution of the Lya absorber 
population conspires to almost cancel out the expected decline in the rms peculiar velocity 
with redshift; hi n = 1 equations (D-(^ yield (Ar2)V2 oc (1 + z)!-^/^ = (1 + 2)0.25_ Qne 
should remember, however, that in difference from galaxies, pressure forces can influence 
the peculiar velocity field of Lya absorption systems, and so naive scalings based on the 
theory of collisionless gravitational instability may not provide accurate predictions for 
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2\l/2 



1.2x10" 



F(0.3,0.7,3) 
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these systems. Since the SZ fluctuations calculated above are dominated by rare high 
column-density systems which are separated by a path length much longer than the 
coherence length of the peculiar velocity field, we are justified in ignoring any velocity 
correlations between neighboring systems on angular scales ~ ^Lya- 
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3. Discussion 

Equation summarizes the expected temperature fluctuations in the microwave sky 
due to the kinematic Sunyaev-Zel'dovich effect of discrete Lya absorption systems. This 
result is based on the observed statistical properties of the Lya forest in the spectra of 
QSOs (cf. Eq. (I)). The predicted anisotropy signal appears on angular scales ~ 0.1-1' and 
decreases roughly as oc on larger scales. Existing VLA observations place an upper 
limit of (AT/T) < 2 x 10"^ on an angular scale of 1.3' (Fomalont et al. 1993), which is 
consistent with the prediction of equation (^. The predicted effect dominates over the 
anisotropics produced by Bremsstrahlung emission from the same Lya absorption systems 
(Loeb 1996) for wavelengths A ^ 3 cm. 

The magnitude of the expected signal from the Lya forest is comparable to the 
O striker- Vishniac effect, which evaluates the anisotropics due to bulk velocities during early 
reionization using second-order perturbation theory (Ostriker & Vishniac 1986; Vishniac 
1987; Hu & White 1995). However, it is important to emphasize that the predictions of 
this Letter rely on the observed population of Lya absorption systems with a minimal 
level of speculation about the full reionization history of the universe. While the thermal 
SZ effect from clusters of galaxies may also add fluctuations of comparable magnitude 
on the same angular scale (Colafrancesco et al. 1994, 1995), it is possible to separate it 
from the kinematic SZ effect of Lya systems either by spectral measurements in several 
microwave bands, or by a cross-correlation analysis with the fluctuations of the X-ray 
sky. The dominant contamination to the signal from Lya systems comes from radio 
sources (see Figs. 9-11 in Tegmark & Efstathiou 1996). Although their contribution to 
the anisotropics on arcminute scales is large, radio sources could in principle be removed 
through high-resolution multiple-frequency observations (Fomalont et al. 1993; Windhorst 
et al. 1993). 

If the anisotropics caused by Lya absorption systems are detected, they would provide 
valuable information about the cosmic velocity fleld and the gas content of these systems 
at high redshifts. The two-dimensional information imprinted on the microwave sky could 
in principle be extended into redshift space, through a cross-correlation analysis with 
absorption spectra of QSOs in the same region of the sky. When extended surveys of 
QSOs, such as the forthcoming Sloan Digital Sky Survey (Gunn & Weinberg 1995), become 
available, it will be possible to perform such a correlation study over large areas of the sky. 

I thank Daniel Eisenstein, Tsafrir Kolatt, and Uros Seljak for useful discussions. This 
work was supported in part by the NASA ATP grant NAG5-3085. 



- 9 - 



REFERENCES 

Bajtlik, S., Duncan, R. C, & Ostriker, J. R 1988, ApJ, 327, 570 
Bechtold, J. 1994, ApJS, 91, 1 

Bechtold, J., Crotts, A. R S., Duncan, R. C, & Fang, Y. 1994, ApJ, 437, L83 

Birkinshaw, M. 1993, in Present and Future of the Cosmic Microwave Background, eds. J. 
L. Sanz, E. Martinez-Gonzales, & L. Cayon, Lecture Notes in Physics, (Springer: 
Berhn), 429, 6 

Blades, J. C, Turnshek, D. A., & Norman, C. A. 1988, QSO Absorption Lines: Probing the 
Universe, (Cambridge: Cambridge University Press) 

Carswell, R. F., Morton, D. C, Smith, M. C, Stockton, A. N., Turnshek, D. A., & 
Weymann, R. J. 1984, ApJ, 278, 486 

Cen, R., Miralda-Escude, J., Ostriker, J. P., & Ranch, M. 1994, ApJ, 437, L9 

Colafrancesco, S., Mazzotta, P. Rephaeh, Y., & Vittorio, N. 1994, ApJ, 433, 454 

. 1995, in the Seventeenth Texas Symposium on Relativistic Astrophysics and 

Cosmology, ed. H. Bohringer, G. E. MorfiU, & J. E. Triimper, Ann. N. Y. Acad. Sci., 
759, 722 

Dinshaw, N., Impey, C. D., Foltz, C. B., Weymann, R. J., & Chaffee, F. H. 1994, ApJ, 437, 
L87 

Dinshaw, N., Foltz, C. B., Impey, C. D., Weymann, R. J., & Morris, S. L. 1995, Nature, 
373, 223 



Fang, Y., Duncan, R. C, Crotts, A. P. S., & Bechtold, J. 1995, preprint |astro-ph/95ToTT2 



Fixsen, D. J., Cheng, E. S., Gales, J. M., Mather, J. C, Shafer, R. A., & Wright, E. L. 
1996, submitted to ApJ, COBE preprint No. 96-11 

Fomalont, E. B., Partridge, R. B., Lowenthal, J. D., & Windhorst, R. A. 1993, ApJ, 404, 8 

Gorski, K., Davis, M., Strauss, M. A., White, S. D. M., & Yahil, A. 1989, ApJ, 344, 1 

Gunn, J. E., & Weinberg, D. H. 1995, in Wide Field Spectroscopy and the Distant Universe, 
ed. S. J. Maddox and A. Aragon-Salamanca (Singapore: World Scientific), in press 



Haehnelt, M. G., & Tegmark, M. 1995, MNRAS, submitted, preprint |astro-ph/9507D77 



- 10 - 



Hernquist, L., Katz, N., Weinberg, D. H., & Miralda-Escude, J. 1996, ApJ, 457, L51 

Hu, W., & White, M. 1995, A& A, submitted, preprint |astro-pli/9507060| 

Janssen, M. A., Scott, D., White, M., Seiffert, M. D., Lawrence, C. R., Gorski, K. M., 
Dragovan, M., Gaier, T., Ganga, K., Gulkis, S., Lange, A. E., Levin, S. M., Lubin, 
P. M., Meinhold, P., Readhead, A. C. S., Richards, P. L., & Ruh, J. E. 1996, ApJ, 
submitted, preprint |astro-ph/960200£ 

Lanzetta, K. M. Wolfe, A. M., & Turnshek, D. A. 1995, ApJ, 440, 435 

Levin, J. J., Freese, K., & Spergel, D. N. 1992, ApJ, 389, 464 

Loeb, A. 1996, ApJ, 459, L5 

Lu, L., Wolfe, A. M., & Turnshek, D. A. 1991, ApJ, 367, 19 

Miralda-Escude, J., Cen, R., Ostriker, J. P., & Ranch, M. 1995, ApJ, submitted, preprint 



Ostriker, J. P., & Vishniac, E. T. 1986, ApJ, 306, L51 

Peebles, P. J. E. 1993, Principles of Physical Cosmology, (Princeton: Princeton Univ. 
Press), 556, 558-559 

Persi, F. M., Spergel, D. N., Cen, R., & Ostriker, J. P. 1995, ApJ, 442, 1 

Pettini, M., Lipman, K., & Hunstead, R. W. 1995, ApJ, 451, 100 

Ranch, M., & Haehnelt, M. G. 1995, MNRAS, 275, L76 

Rephaeh, Y. 1995, ARA&A, 33, 541 

Sargent, W. L. W., Young, P. J., Boksenberg, A., & Tytler, D. 1980, ApJS, 42, 41 

Sargent, W. L. W., Steidel, C. C, & Boksenberg, A. 1989, ApJS, 69, 703 

Sembach, K. R., Steidel, C. C, Macke, R. J., & Meyer, D. M. 1995, ApJ, 445, L27 

Smette, A., Robertson, J. G., Shaver, P. A., Reimers, D., Wisotzki, L., & Kohler, Th. 1995, 
A&A Suppl., 113, 199 

Spitzer, L. 1978, Physical Processes in the Interstellar Medium, (New- York: Wiley), 107 

Stengler-Larrea, E. A., Boksenberg, A., Steidel, C, Sargent, W. L. W., Bahcall, J. N., 
Bergeron, J., Hartig, G. P., Jannuzi, B. T., Kirhakos, S., Savage, B. D., Schneider, 
D. P., Turnshek, D. A., & Weymann, R. J. 1995, ApJ, 444, 64 



- 11 - 



Storrie-Lombardi, L. J., McMahon, R. G., Irwin, M. J., & Hazard, C. 1994, ApJ, L13 

Strauss, M., Cen, R., & Ostriker, J. R 1993, 408, 389 

Strauss, M. A., & Willick, J. A. 1995, Physics Reports, 261, 271 

Sunyaev, R. A., & Zel'dovich. Ya., B. 1980, ARA&A, 18, 537 

Tegmark, M., & Efstathiou, G. 1996, MNRAS, in press, preprint |astro-ph/9507009 
Tytler, D. 1987, ApJ, 321, 49 

Walker, T. R, Steigman, G., Kang, H., Schramm, D. N., & Ohve, K. A. 1991, ApJ, 376, 51 
Webb, J. K., & Barcons, X. 1991, MNRAS, 250, 270 

Windhorst, R. A., Fomalont, E. B., Partridge, R. B., & Lowenthal, J. D. 1993, ApJ, 405, 
498 

Wright, E. L., Mather, J. C., Fixsen, D. J., Kogut, A., Shafer, R. A., Bennett, C. L., 
Boggess, N. W., Cheng, E. S., Silverberg, R. P., Smoot, G. P., & Weiss, R. 1994, 
ApJ, 420, 450 

Zhang, Y., Anninos, P., & Norman, M. L. 1995, ApJ, 453, L57 



This preprint was prepared with the AAS M5t;X macros v3.0. 



